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ABSTRACT 

We present new information on galaxies in the vicinity of luminous radio galaxies and quasars at 
z ~4, 5, and 6. These fields were previously found to contain overdensities of Lyman Break Galaxies 
(LBGs) or spectroscopic Lya emitters, which were interpreted as evidence for clusters-in-formation 
( "protoclusters" ) . We use HST and Spitzer data to infer stellar masses from stellar synthesis models 
calibrated against the Millennium Run simulations, and contrast our results with large samples of 
LBGs in more average environments as probed by the Great Observatories Origins Deep Survey 
(GOODS). The following results were obtained. First, LBGs in both overdense regions and in the 
field at z = 4—5 lie on a very similar sequence in a z'— [3.6] versus 3.6/im color-magnitude diagram. This 
is interpreted as a sequence in stellar mass (M* ~ 10 9 -10 n Mq) in which galaxies become increasingly 
red due to dust and age as their star formation rate (SFR) increases, while their specific SFR stays 
constant. Second, the two radio galaxies are among the most massive objects (M* ~ 10 11 Mq) known 
to exist at z ~ 4 — 5, and are extremely rare based on the low number density of such objects as 
estimated from the ~25x larger area GOODS survey. We suggest that the presence of the massive 
(radio) galaxies and associated supermassive black holes has been boosted through rapid accretion of 
gas or merging inside overdense regions. Third, the total stellar mass found in the z = 4 protocluster 
TN1338 accounts for <30% of the stellar mass on the cluster red sequence expected to have formed 
at z > 4, based on a comparison with the massive X-ray cluster CU252 at z = 1.2. Although future 
near-infrared observations should determine whether any massive galaxies are currently being missed 
by our UV/Lya selections, one possible explanation for this mass difference is that TN1338 evolves 
into a smaller cluster than C11252. This raises the interesting question of whether the most massive 
protocluster regions at z > 4 remain yet to be discovered. 

Subject headings: galaxies: high-redshift — galaxies: stellar content — (cosmology:) large-scale struc- 
ture of universe — galaxies: individual (TN J1338-1942,TN J0924-2201) — (galax- 
ies:) quasars: individual (SDSSJ0836+0054,SDSSJ1030+0524) 



1. INTRODUCTION 

Galaxy clusters formed in high-density regions of dark 
matter that coll apsed earlier than surrounding low- 
density regions (|Kaiserlfl98l ISpringel et al.ll2005l ). and 
models predict that the galaxies that formed inside these 
high-density regions also may have formed earlier or 
evol ved at an increased rate compared to their surround- 
ings (|Kauffmannlll995tlBenson et alj2001l:lThomas et~all 
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120051 : iDe Lucia etaD |2006). The mass functions of 
clusters at low and high redshift can significantly im- 
prove current estimat es of the cosmological parameters 
(|Vikhlinin et al.ll2009f ). Clusters of galaxies are thus im- 
portant not only for the formation of dark matter struc- 
ture, but also for studying processes related to galaxy 
formation, such as the history of gas accretion and deple- 
tion, the feedback of supernovae and active galactic nu- 
clei (AGN), the enrichment of the intracluster medium, 
and the spectra-morphological transitions of galaxies. 

In the local Universe, the galaxies in the centers of 
clusters inhibit a distinct location in the color-magnitude 
diagram (CMD). This cluster "red sequence" consists of 
predominantly spheroidal and lenticular galaxies having 
old stellar populations and high stellar masses, and usu- 
ally includes the brightest cluster galaxy. In order to 
understand how and when this relation formed, it is im- 
portant to study clusters at high redshifts. Traditional 
surveys have uncovered clusters at redshifts as high as 
z ~ 1.5 either by searching for galaxies on the character- 
istic red sequence or by searching for hot cluster gas in 
the X-rays. A red sequence of morphological early-type 
galaxies is prominently present in even the most distant 
clusters, and the small scatter in galaxy colors around 
the sequence has been used to infer formation redshifts 
of ~ 2-5 (e.g lEllis et al.lll997tlva"n Dokkum et~ai] | 2000 | : 
iStanford et "all 120021 120051 : iBlakeslee et all 12003d 120061 : 
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Holden et all 120051: i Postman et all 120051 : iMullis et all 
2005 HMei et alJl200a iRettura et alll2008| ). 
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The epoch corresponding to the assembly of massive 
cluster galaxies is presumed to be marked by a violent 
stellar mass build-up. In order to target this epoch, var- 
ious techniques have been used. One technique is based 
on the hypothesis that at very early times the most mas- 
sive galaxies trace the highest density regions and can 
thus be used as "lighthouses" to pinpoint possible "pro- 
toclusters" . Radio galaxies (RGs) are among the largest 
and mos t massive galaxies at z = 2 — 5 (M* ~ 10 10 ~ 12 
M R , e.g. Pcntericci et al. 1997; Rocca-Volmcrange et al. 
12004 ISevmour et~afl 120071 : lHatch et al.l I2009L this pa- 
per), while the most luminous high rcdshift quasars 
(QSOs) have some of the largest black hole masses 
inferred (M RH ~ lO 8 " 10 M Q , e.g. Uiang et all [20071 : 
iKurk et al1l2007h . Overdensities (of the order of a few) 
of star-forming galaxies suggest that at least some ra- 
dio galaxies are asso ciated with rare peaks i n the large- 
scale structure (e.g. | Pentericci et al.l 120001 : iMilev et alJ 
12004 IVenemans et all I2002L 120071 ). Similarly, the num- 
ber counts of faint z-dropouts towards two z ~ 6 QSOs 
(jStiavelli et al.ll2005l : iZheng et al]l200?i iKim et al.ll2009h 
also i ndicate a possible enhancement (see lOverzier et all 
2009, for a detailed discussion). Similar overdensities 
have also been found in "rand o m" fields at z ~ 2 — 6 
(e.g. ISteidel et all I1998L 120051: iShimasaku et al] 120031 : 
lOuchi et alJ 1200a lOta et al.ll2008l ). and are consistent 
w ith being pro g enitor s of present-day clusters. 

ISteidel et all (pOOfih found an increase in the stel- 
lar masses and ages for star-forming galaxies in a pro- 
tocluster at z = 2 . 3 com pared to the surrounding 
field. InE irm et all (|2008l ) we presented evidence for 
a population of relatively evolve d galaxies near a ra - 
dio galaxy at z = 2 2 (see also iKodama et al.l I2007D . 
and IKurk et al.1 (|2004 ) found evidence for different spa- 
tial distribution s of b lue and red galaxies. Although 
IVenemans et al.l ([2007D found no differences in the phys- 
ical properties of LAEs in protoclusters compared to the 
field, prot oclusters tend to have an elevated AGN frac- 
tion (e.g. IPentericci et al.1 120021: IVenemans et al.l 120071 : 
ILehmer et al.ll2008D . iLehmer et al.l ([20081) found an en- 
hancement in the stellar ma sses of LBGs in t he pro- 
tocluster SSA22 a t z = 3.1 (ISteidel et al.l I1998T) . while 
iPeter et al.1 (|2007h and lOverzier et alJ (|2008l ) did not 
find any morphological differences between UV-selected 
galaxies in protoclusters compared to the field using data 
from the Advanced Camera for Surveys (ACS) on-board 
the Hubble Space Telescope (HST). 

The ACS, in particular, has enabled the selection of 
lar ge samples of fa int UV "dropout" gal axies at z > 
4 (|Giavalisco et al.l 12004 iBouwens etabl 12004 I2007D . 
Combined with the great sensitivity in the mid-infrared 
provided by the Infrared Array Camera (IRAC) on 
the Spitzer Space Telescope (SST), w e can study the 
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stellar population s of those galaxie s (lYan et al 
Labbe et al. l 120061: lEvles et al.l 120071: IStark et al.l 12007 
20091 : IVerma et al.ll2007l : lYabe et al.ll2008l ). While obser 



vations in the UV probe the (dust attenuated) light from 
young, massive stars, the rest-frame optical light offers a 
more direct probe of the total stellar mass. Results at z ~ 
2 — 3 imply that some LBGs and near-infrared selected 
galaxies were alread y quite massive at z > 5 implying 
rapid evolution (e.g. IShaplev" et alJl2001t iPapovich etall 



20071 ). Recently. IStark et al.l (|2009h found that the z > 4 
UV luminosity function (LF) is dominated by episodic 
star formation (<500 Myr on average) from systems 
that reached their high UV luminosities in short times 
<300 Myr), consiste nt with the decline in the UV LF 
Bouwens et al.ll2007l ). Galaxies having relatively high 
equivalent widths of Lya emission (the so-called "Lya- 
emitters") have received attention (partly due to the 
relative ease of their detections out to the highest red- 
shifts) concerning the question whether they mark a par- 
ticular phase in the evolution of typica l starbursts (e.g. 
I Shaplev et"aLll2003l:lGawiser et alJ l2006b IPente ricci et all 
120071 12009t lOuchi et al.ll2008T) . 

In this paper we will derive stellar masses, as inferred 
from the rest-frame optical fluxes, for a sample of LBGs, 
LA Es and RGs 9 in two protoclusters: "TN1338" at z ~ 
4.1 (IVenemans et al.|l2002t ;lMilev et~aLll2004 IZirm et all 
120051: lOverzier et al.1 12001. and "TN0924" at z ~ 5.2 
(jVenemans et al.l [2004 lOverzier et al.l [20061 ) . For com- 
pleteness, we also present a small sample of i-dropouts 
near t wo q uasars at z ~ 6 s tudied by iStiavelli et al.l 
(|2005h and IZheng et al.l (|2006f ). The target fields are 
listed in Table [TJ Comparison with large samples from 
the GOODS survey (jGiavalisco et al.ll2004l) will enable us 
to compare the overdense regions with more average en- 
vironments. While the data presented here cannot com- 
pete in terms of eit her depth or area with other studies in 
the literature (e.g. IStark et al.ll20 09). our data is unique 
in the sense that we cover a few sightlines towards rare, 
overdense regions and luminous AGN. Specifically, we 
will address the following questions: 

(1) How do the stellar masses of radio galaxies, LBGs 
and LAEs compare? 

(2) Do galaxies in overdense regions at high redshift have 
distinct physical properties? 

(3) Is the total amount of stellar mass detected in pro- 
toclusters consistent with the stellar mass in galaxies on 
the low redshift cluster red sequence? 

(4) Have we missed a population of obscured or quiescent 
galaxies not selected by our UV/Lya selections? 

The structure of this paper is as follows: The samples 
and our measurement techniques are described in §2. In 
§3 we present the rest-frame UV-optical color-magnitude 
diagrams at z ~ 4, 5,6 and estimate stellar mass distri- 
butions based on simple model comparisons. We dis- 
cuss our results in §4. All magnitudes are in the AB 
system. For comparison with other studies, it is conve- 
nient to s pecify that the char acteristic luminosity L* =3 
at z ~ 3 (|Steidel et al.l I1999L i.e., M 1700:AB = -21.07) 
corresponds to z-band magnitudes of ^24.9, 25.3 and 
25.6 at z^4, 5 and 6, respectively. We adopt the fol- 
lowing cosmological parameters: J7m = 0.27, J7a = 0.73, 
and H = 73 km s _1 Mpc -1 . 



2. SAMPLES, DATA AND MEASUREMENTS 
2.1. Primary Samples 
2.1.1. Radio galaxy fields 

9 Stellar mass estimates for the r adio galaxies TNJ133 8-1942 
and TNJ0924-2201 are also given in ISevmour et al.l d2007h . Our 
analysis differs in that it is based on our deeper data and uses more 
sophisticated photometry. However, our results are consistent. 
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This paper is based on various samples of galaxies in 
overdense regions selected as part of previous works (see 
§1). For the two radio galaxy fields we make use of large 
samples of LBGs and LAEs selected from our previous 
HST/ACS imaging and Lya surveys. 

In the field TN1338 we will focus on a sample of LAEs 
and LBGs selected near the radio galaxy TN J1338-1942 
at z = 4.11. The field covered wi th both HST and Spitze r 
contains 12 Lya emitters from [Vcncma ns~et al.l (2007): 
these LAEs are part of a larger structure that is « 5 ± 1 
times more dense compared to average expectations, and 
has a velocity dispersion of w 300 km s^ 1 centered on the 
radio galaxy. The LBG sample consists of 66 ^475-band 
(g ) dropouts brighter t han z 8 50=27.0 mag. As detailed 
in lOverzier et a l. (2008), the photometric redshift distri- 
bution peaks at z = 4.1, and has a width of Az ~ 0.5 
(FWHM). About half of the LAEs (\z - Az\ < 0.03) 
are also selected as part of the dropout sample, indicat- 
ing that the true width of the redshift distribution of 
the dropouts could be narrower than suggested by the 
photo- z analysis given its fairly limited accuracy. Unfor- 
tunately, the redshifts of most of the dropouts are diffi- 
cult to determine spectroscopically due to their faintness 
and because all high equivalent width Lya sources have 
already been identified. Statistical studies indicate that 
the ratio of Lya-faint to Lya-bright LBG s in a UV flux - 
limited sample is about 3:1 (e.g. iShapley et ai1l2003f ). 
The population of g-dropouts represents an overdensity 
at the > 5a significance with respect to the average field 
as estimated from the GOODS survey, consistent with a 
substantial fraction lying near the redsh ift of the radio 
galax y and its associated Lya emitters ()Overzier et al.l 
2008). The LBG/LAE samples are listed in Tables M 
and El 

In the field TN0924 we make use of our samples of 
LAEs and LBGs selected near the radio galaxy TN 
J0924-2201 at z = 5.20. This galaxy is the highes t red- 
shift radio galaxy known (jvan Breugel et ali ri999L The 
radio galaxy has six spectrosco pically confirmed LAE 
companions (\z — Az\ < 0.015: IVenemans et alJ f2004). 
four of which are in the field covered by HST and Spitzer. 
The LBG sample consists of 23 V(306-band (V) dropouts 
brighter than Z85o=26.5 mag. Analogous to the TN1338 
sample, the U-dropout selection was tuned to the se- 
lection of candidate LBGs n ear the redshift of the radio 
galaxy (|Overzier et al.ll2006f) . The TN0924 field contains 
about twice as many U-dropouts compared to the aver- 
age number found in GOODS, with a ~1% chance of 
finding a similarly high number also based on a compar- 
ison with GOODS. Their photometric redshift distribu- 
tion has a width of Az ~ 0.5 (FWHM), but its true width 
could be narrower based on the same arguments as given 
above. Two of the LAEs are also selected as dropouts, 
implying that, out of all candidates selected, the actual 
number of LBGs physically associated with the radio 
galaxy and other LAEs may not be higher than ~ 8 for a 
Lya faint-to-bright ratio of 3:1 (see lOverzier et aT1l2006l . 
for discussion) . The samples are listed in Table |U 

2.1.2. Quasar Fields 

For the QSO field SDSS0836 we make use of a sample of 
j-dropouts that was selected in the 3.4'x3.4' HST/ACS 
pointing towards the radio-loud QSO SDSS J0836+0054 
at z = 5.82. This QSO is the highest redshift radio- 



loud quasar known. I Zheng et al.l (2006) found 7 dropouts 
within a 5 arcmin 2 projected region near the QSO, hav- 
ing colors consistent with them being at z sa 5.8. In ad- 
dition, we obtained archival ACS observations (G09777, 
PI: Stiavelli) of the field towards QSO SDSS J1030+0524 
at z = 6.28 {QSO field SDSS10 30). The data was pro - 
cessed with the pipeline APSIS (jBlakeslee et al.ll2003bf) . 
and we perf ormed our own selection of i-dropouts using 
SExtractor (iBertin fc Arnoutsll996l) and using the same 
criteria as in lZheng et al.l ( 20061 ): i — z > 1.3 (magJso), 
z < 26.5 mag, and a SExtractor star/galaxy index < 0.8. 
We visually inspected the ACS images and excluded can- 
didates near the detector edges. Because of our own 
processing of the data as well as more stringent selec- 
tion criteria, the number of our candidates (5) in the 
SDSS 1030 field is somew hat smaller than tha t given pre- 
vious ly in the literature (|Stiavelli et al.ll2~005l : iKim et all 
2009). The two samples of i-dropouts in QSO fields are 
listed in Table El 

2.2. GOODS Comparison Data 

Our comp arison data i s base d on the large dropou t 
samples from lStark et alj l|2007l) and lStark et al] ([2009), 
who performed a detailed study of the high redshift 
LBG populations in the ~ 320 arcmin 2 GOODS sur- 
vey (|Giavalisco et al.l 12004). The parent samples from 
Stark et al. consist of 2443, 506 and 137 B-, V-, and i- 
dropouts to a limiting optical magnitude of 28, of which 
^35% was found to be sufficiently isolated in the IRAC 
images for accurate photometry. The selection criteria 
for these samples are equivalent to those used for our 
protocluster fields. Details on the selection and the HST 
and Spitzer phot ometry of the GOODS samples can be 
found elsewhere (jStark et al.ll2007l I2009D . 

For comparison with our TN1338 field, we will make 
use of their large sample of z w 4 _B-dropouts, while for 
TN0924 we will use the sample of z » 5 U-dropouts. 
For the two QSO fields at z w 6 we use a small sample 
of i-dropouts also from GOODS. After applying appro- 
priate magnitude cuts, we are left with 420, 122 and 42 
comparison objects, respectively. Because the overdense 
fields listed above have all been previously demonstrated 
to be overdense relative to the much larger GOODS sur- 
vey, the latter is expected to provide a good handle on 
the LBG population in more average environments. 

2.3. Millennium Run Simulations Comparison Samples 

In order to test our strategy for determining stellar 
masses based on a comparison between simple SED mod- 
els and the observed magnitudes, we have selected a large 
sample of mock galaxies simulated as p art of the Mil- 
lenni um Run (MR) Simulations project (Springel et al. 
|2005[ ) . In brief, galaxies were modeled by applying semi- 
analytic model (SAM) prescriptions of galaxy formation 
to the dark matter halo merger trees from the MR. The 
SAM includes processes such as gas cooling, star for- 
mation, reionization heating, and feedback from super- 
nova and AGN that are coupled to the stellar synthesis 
model library from BC03 in order to calculate the re- 
sulting spectrum of each galaxy in the simulations as its 
evolution is tracked across discrete "snapshots" in time 
between z = 127 and z = (see ICroton et al.l 120061 : 
iDe Lucia fc Blaizotl [20071 : IGuo fc White! 120081 . for de- 
tails). We use the mock galaxy catalogues 'delucia2006a' 
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Fig. 1. — Panels show the ACS zsso-band and IRAC 3.6/im 
images of the radio galaxies TN1338 at z = 4.1 (bottom panels) and 
TN0924 at z = 5.2 (top panels). The images measure 20"x20". 
The ACS images were convolved with a Gaussian kernel of O'.'l 
(FWHM). The IRAC photometry was obtained by deblending the 
radio galaxies and their neighbouring sources using GALFIT (see 
§2 for details). 

from iDe Lucia fc BlaizotJ ((2007) and randomly select 1 
million galaxies each at z « 4 and z w 5. Previous stud- 
ies have shown that the simulations are able to reproduce 
the observed p roperties of galaxies at these redshifts rela - 
tively well fe.g. lGuo fc Whitdl2008l : [Overzier et al1l2009h . 
The large size of the mock galaxy catalogues selected en- 
sures that we sample a wide range in SFHs, ages, masses, 
dust attenuations, metallicities and large-scale environ- 
ments. 

2.4. IRAC Imaging 

Follow-up infrared images for the fields TN0924 and 
SDSS0836 were obtained with the IRAC on the SST 
(PID 20749, Zheng 2006). We retrieved archival IRAC 
data of TN1338 (PID 17, Fazio 2004) and SDSS1030 
(PID 30873, Labbe 2007). The IRAC camera includes 
imaging data at 3.6, 4.5, 5.8, and 8.0 /im (channels 1-4), 
each with a 5.'2 x 5f 2 field and a pixel size of ~ 1'.'22. 
We concentrate here on the shortest-wavelength Spitzer 
image (channels 1 and 2) for the field at z ~ 4 (TN1338), 
and on the 3.6/rni images for the higher redshift fields. 
The exposure times at 3.6/^m are 17 ksec (TN0924 & 
SDSS0836) and 10 ksec (TN1338 & SDSS1030). For our 
photometry, we used the pipeline-processed IRAC im- 
ages at the 'post-basic calibrated data' (PBCD) stage. 
The main steps in the pipeline include dark current sub- 
traction, flat fielding, flux calibration (in units of MJy 
sr _1 ), geometric distortion correction and mosaicking of 
the individual frames. Sensitivities of each image were 
estimated from the 2 a standard deviation in the flux 
distribution measured in randomly placed 3" diameter 
apertures. The 2a limiting AB magnitudes are 25.5 
(TN0924), 25.2 (J0836), 25.2 (TN1338), and 25.2 (J1030) 
mag in the 3.6 /im channel, and 25.0 (TN1338) and 24.6 
(J1030) mag in the 4.5 /im channel. 



2.5. IRAC Source Photometry using GALFIT 

The IRAC images suffer from overcrowding due to its 
relatively large point spread function (PSF) of ~ l'.'S 
(FWHM) and great sensitivity. To address this is- 
sue, we use a deblending technique whereby contaminat- 
ing n eighbors are subtracted using GALFIT (|Peng et all 
2002) by performing a fit to the objects of our interest 
and all their close neighbours simultaneously. GALFIT 
constructs a two-dimensional model of the data according 
to specified input parameters (e.g., positions, magnitude, 
effective radius, axis ratio, and position angle), performs 
a convolution with the instrument PSF, and fits the re- 
sult to the data through an iterative ^-minimization 
process. The IRAC PSF is obtained directly from the 
images by stacking 10 bright, isolated point sources. We 
use the higher resolution ACS z-band image as the refer- 
ence for the initial GALFIT input parameters. GALFIT 
was then applied to the IRAC images. Two versions of 
the fitting process were carried out. First, all input pa- 
rameters as determined from the z-band image were held 
fixed, except the magnitudes. The resulting \ 2 value 
and residual image was then visually inspected to deter- 
mine the goodness of fit. When the best possible fit was 
not achieved, we repeated the GALFIT procedure, but 
this time allowing all input parameters to vary. Each 
object was individually and interactively processed until 
the most successful fit was achieved, while avoiding the 
oversubtraction of the flux contribution from the neigh- 
boring sources. Following this, the two-dimensional mod- 
els of the light profiles of the LBGs/LAEs and nearby 
contaminating sources were deblended. This deblend- 
ing process works quite well for determining the intrinsic 
fluxes of most of LBGs/LAEs. Those sources for which 
GALFIT failed to satisfactorily deblend the emission of 
LBGs/LAEs from the neighbour sources were removed 
from the sample. This leaves a total of 60 sources (38 
LBGs and three LAEs at z ~ 4.1, 15 LBGs at z ~ 5.2 
and 4 LBGs at z ~ 6), for which reliable photometry 
from GALFIT was obtainable. 

The measured 3.6 (4.5) /jm magnitudes for the four 
fields are listed in Tables 1-4, along with the classifica- 
tions of galaxies for their detections and degree of con- 
fusion, as (1) isolated and detected; (2) isolated but un- 
detected; (3) confused (4) heavily blended. In Fig. [T] 
we show the ACS zsso and IRAC 3.6/im images of the 
two bright radio galaxies for an illustration of the data 
quality and blending issues. 

2.6. Emission line contamination corrections 

Because the IRAC 3.6/im flux of the radio galaxy TN 
J1338-1942 is significantly contaminated by Ha flux from 
its large emission line halo that is typical of RGs, we have 
used the emission-line free flux at 4.5/im and a small 
model extrapolation based on the typical [3.6]-[4.5] col- 
ors of star-forming galaxies at z « 4 to estimate the true 
continuum magnitude at 3.6/im as follows. The mea- 
sured [3.6]-[4.5] color was -0.65 mag (see Table 1), while 
a comparison with stellar population models (see below) 
in 4.5/mi vs. z-[4.5] space would predict a color of [3.6]- 
[4.5]~0-0.2 mag. We thus assume that Ha contributes 
about 0.75 ± 0.10 magnitude to the 3.6/im measurement 
of the radio galaxy, and we will apply the correction to 
all subsequent figures and tables. The 3.6/im flux of the 
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Fig. 2. — Transmission curves of the relevant optical (HST/ACS), 
near-infrared, and infrared (IRAC) filters. We overplot several 
example spectral energy distributions modeled using BC03. Lines 
show the expectation for an exponentially declining star formation 
history with r = 100 Myr at ages of 10 (dashed lines) and 100 
Myr (solid lines), and redshifted to z = 4.1 (blue lines), z = 5.2 
(green lines), and z = 6.0 (red lines). A small attenuation of 
E(B — V) = .1 magnitude was modeled using the prescription of 
Calzctti (2001) and applied to each of the tracks. 

radio galaxy TN J0924-2201 is not affected by Ha given 
its higher redshift of z = 5.2. 

While the 3.6/j,m flux of z ~ 4 star-forming galaxies 
such as LBGs and LAEs is expected to have an enhance- 
ment due to Ha as well, we do not correct for this effect 
as we can assume that galaxies in both GOODS and in 
the protocluster fi elds will be affected in a similar way. 
IStark et al.l (|2007h estimate that Ha could in principle 
increase the broad band flux by ~ 10 — 20%. 

2.7. Estimating stellar masses 
2.7.1. Rationale 

In order to answer the question of whether there 
are any differences between the different types of 
galaxies across the CMD and across environment, 
we will seek to determine estimates of (primarily) 
stellar masses based on a set of z— [3.6] versus 3.6/im 
CMDs. Only recently has it become feasible to attempt 
an interpretation of the observed magnitudes of the 
relatively faint galaxy populations at z ~ 3 — 6 in 
terms of stellar masses, ages, dust and star formation 
histories based on detailed fitting of SEDs using model 
libraries (see references given in §1 and therein). While 
there is ongoing debate about the correctness of the 
treatment of different stell ar evolutionary phases in 
the model libraries (e.g. [Fioc fe Rocca-Volmerangel 
1997b iBruzual fe Chariot] 120031; iMaraston et all 120061: 
Rettura et al.l 120061 : iKannappan fe Gawiserl 120071 : 
Eminian et al.l l2008f) . they at least allow for an as- 
sessment of the comparative properties among samples 
across redshift or across the CMD, provided that these 
samples are sufficiently 'simple' to model and that the 
adopting of a faulty model does not lead to catastrophic 
systematic effects across the various samples. At the 
least, this requires highly accurate photometry in all 
of the rest-UV, optical and infrared, and such studies 
have therefore been limited to a small set of well-studied 
fields. 

While our ACS imaging is of comparable depth to, 
e.g., GOODS, we lack sufficiently deep data in the ob- 
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Fig. 3. — The i— z versus z-[3.6] color-color diagram for the 
TN1338 field. Symbols indicate the radio galaxy (large filled cir- 
cle), LBGs (open circles) and LAEs (filled stars). Tracks indicate 
the colors of models having r = 30 Myr (dotted lines) and r = 200 
Myr (solid lines) with ages in Myr indicated along the tracks. At- 
tenuations of E(B — V) = 0.0,0.15,0.4 mag are shown in blue, 
green and red, respectively. Most objects do not require a large 
amount of reddening due to dust. 

served near-infrared. Our IRAC imaging provides a rela- 
tively good match to the ACS imaging (most UV-selected 
sources detected), but it is not as deep as some of the 
large field surveys. Therefore, we do not believe that a 
detailed SED fitting is warranted by our data, nor is it 
necessary for the limited questions we wish to ask. In- 
stead, as we will show below, we can obtain sufficiently 
accurate information on stellar masses based on the com- 
parison with simple baseline models rather than SED fit- 
ting ( £|2.7.2[) . We will test our method by making a com- 
parison with the properties of mock galaxies selected in 
a large cosmological simulation ( q2.7.3[) . 

2.7.2. Comparison with Stellar Synthesis Models 

As detailed in lStark etail (|2009h . at z > 4 the 3.6/xm 
channel still probes only as red as the rest-frame opti- 
cal, where the contribution from TP-AGB stars is much 
less severe than in the near-infrared. As a result, the stel- 
lar m ass estimates based on "BC03" (|Bruzual fe Chariot! 
2003; models are expected to be higher by at most ~30% 
when compared to the newer gen eration of models that 
include this phase (jBruzua l [2007Tl . and the effect of AGB 
stars becomes less severe with increasing redshift and 
metallicity with respect to subsolar metallicity models 
at z ~ 4. Furthermore, the SED fitting appears to be 
relatively insensitive to whether redshifts are being fixed 
or allowed to vary within their Az — ±0.25 uncertainty 
distributions typical for dropout samples. 

In order to interpret the data in the CMDs that will 
be presented in Sj31 we have selected a set of BC03 mod- 
els that we will use for a simple baseline comparison 
with our data. The models are based on the best-fit 
ages as obtained from pan-chromatic SED fits performed 
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bv lStark et al.l (|2009), and the typical atte nuations due 
to dust as found bv iBouwens et alJ (|2009[ ). Both stud- 
ies are based on large samples of dropout galaxies at 
4 ^$ z < 6 for which deep photometry is available. We 
choose an exponentially declining star formation history 
with an e-folding timescale r of 100 Myr. We consider a 
range of ages from 100 to 300 Myr, straddling the best- 
fit age of -200 Myr as found bv lStark et all |2009D . We 
will also consider three levels of attenuation by dust in 
terms of E(B — V) =0.0 mag (no dust), 0.15 mag (mod- 
erate extinction) and 0.3 magnitude (dusty), modeled 
using the recipes for starburst galaxies given in lCalzettil 
(2001). All models have a metallicity similar to that of 
the Large Magellanic Cloud. We convert the models to 
the appropriate redshift, and calculate colors and mag- 
nitudes by convolving the spectral energy distributions 
with the HST and Spitzer filter transmission curves (Fig. 
[2]). Given an age and an attenuation, an estimate for the 
stellar mass of each object is obtained by determining 
the absolute mass scale of the model required in order to 
reproduce the 3.6/ im magnitude observe d. Motivated by 
a recent study by IBouwens et al.l (|2009t ) that finds clear 
evidence of correlations between the UV color and abso- 
lute UV luminosity of B- and V— dropout galaxies, we 
will also investigate a model in which the dust and age 
are allowed to vary as a function of absolute magnitude. 

These simple model comparisons show that while the 
3.6/im magnitude is an approximate gauge of stellar mass 
at each of the redshifts, the z-[3.6] color is affected by 
significant degeneracies between age, dust and SFH that 
can only be satisfactorily addressed using deeper data 
and a wider set of bands than currently available for our 
fields. However, since this paper mainly focuses on gen- 
eral trends in stellar mass for different populations, our 
main results will not be significantly affected by this is- 
sue. Nonetheless, for the TN1338 field we can be a little 
more specific and try to get a better handle on the pres- 
ence of dust. The i— z UV continuum color reacts more 
strongly to reddening due to dust than the z-[3.6] color. 
This is illustrated in Fig. [3] where we show the TN1338 
i—z vs. z~[3.6] color-color diagram compared to mod- 
els with three different attenuations of 0.0 (blue lines), 
0.15 (green lines) and 0.4 mag (red lines). The sym- 
bols indicate LBGs (open circles), LAEs (filled stars), 
and the RG (filled circle). Most objects are very blue 
in the UV (i-z « 0.0 mag) implying a small amount of 
attenuation at most. In other words, strongly reddened 
galaxies do not appear to make an important contribu- 
tion. This is not surprising since not only is our UV 
selection biased against very dusty objects, recent evi- 
dence suggests th at dusty objec ts become increa singly 
rare at z > 4 (e.g. iMcLure et al.l l200l IStark et al.1120091 ; 
IBouwens et al.ll2009l Sect. [4T4lof this paper). 

2.7.3. Comparison with Mock Galaxies in the Millennium 
Run Simulations 

Above we have outlined our basic strategy of using sim- 
ple baseline models to convert 3.6/zm fluxes into stellar 
mass estimates. In order to test this we use our large 
sample of mock galaxies from the Millennium Run sim- 
ulations (see Sect. 12. 3[) . In the top-left panel of Fig. [4] 
(panel a) we plot the stellar mass against the observed- 
frame 3.6/xm magnitude for each of the mock galaxies at 
z = 4. The 3.6/im magnitude is clearly correlated with 
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Fig. 4. — The accuracy of our method of retrieving stellar masses 
based on simple SED assumptions, tested against one million mock 
galaxies selected from the Millennium Run (MR) cosmological sim- 
ulations. Top panels show the relation between stellar masses and 
3.6/im magnitudes for the entire MR sample (panel a, left) and 
for a cleaned subset of galaxies most similar to our UV-selected 
samples (panel b, right). Blue points indicate the median val- 
ues and the lcr standard deviation in bins of Am = 0.3 mag. A 
second-order polynomial fit to the correlation shown in panel b 
is logM*[M ] = 11.51 + 0.32m 3 . 6 - 0.016m§ 6 (red line). In the 
bottom panels, we demonstrate how accurately we can derive the 
true stellar masses from the observed 3.6/im magnitudes under the 
assumption of a fixed baseline SED characterized by r = 100 Myr, 
t = 200 Myr, E(B - V) = 0.15 mag, and LMC metallicity. Despite 
the inevitable mismatches between the baseline model and the true 
properties of the simulated galaxies, there is still a very good cor- 
relation between the true mass and the mass inferred using the 
baseline model. The one-to-one correlation is indicated by the red 
dashed lines. There is a small systematic offset in the sense that 
for more massive galaxies we underestimate their true mass based 
on the baseline SED model (by about ~0.3 dex at M, ~ 10 11 Mq, 
see panel d for the sample most similar to our UV selection). 



the stellar mass, indicating that the flux at rest-frame 

~7000A is a good tracer of the mass. Next, we limit 

the mock sample to only those galaxies that are most 

similar to the galaxies in our UV-selected LBG samples 

by requiring relatively little dust (E(B — V) < 0.3 mag, 

see Fig. [3]) and current star formation activity at a rate 

of >1 Mq yr _1 . The result is shown in the top-right 

panel (b). The scatter around the correlation between 

M* and 3.6/im is significantly reduced. We see that a 

range in the 3.6/im magnitude from 26 to 21 mag cor- 

li 



0.016m| 6 (red line in 



responds to a range in stellar mass from 10 a to 10 
Mq. A second-order polynomial fit to the data yields 
logM*[M Q ] = 11.51 + 0.32m 3 . 6 
panel b) 

As stated earlier there are ambiguities between SFH, 
age and dust that cannot be solved for based on our cur- 
rent data. However, we can also use the MR simulations 
to test how well we can retrieve stellar masses when as- 
suming that each galaxy can be approximated by a single 
one of our simple baseline SEDs from Sect. 12.7.21 (e.g., 
r = 100 Myr, t = 200 Myr, E(B - V) = 0.15 mag) for 
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each 3.6/im magnitude. The result is shown in the bot- 
tom panels of Fig. |4j again for all (left panel, c) and for 
the "cleaned" subsample (right panel, d). Despite the in- 
evitable mismatches between the baseline model and the 
true properties of the simulated galaxies, there is a very 
good correlation between the true mass and the mass 
inferred using the baseline model. We do however see 
a systematic offset in the sense that for more massive 
galaxies we underestimate their true mass based on the 
baseline SED model (by about -0.3 dex at M* - 10 11 
Mq). This can be understood because more massive ob- 
jects have higher SFRs and have more dust in the sim- 
ulations. Therefore, they also will have more dust com- 
pared to our baseline model so their 3.6/im flux is lower 
and hence we get a lower mass using the baseline model 
conversion. 

The results obtained for our z = 5 sample are very 
similar, but are not shown here for brevity. In summary, 
these simulations demonstrate that our approach of using 
simple baseline models to infer stellar masses from the 
3.6/im data is justified. 

3. RESULTS 

In this Section, we present the 3.6/im versus z— 
[3.6], or rest-frame UV-optical, color-magnitude dia- 
grams (CMDs) of the various galaxy populations de- 
tected in the RG and QSO fields. We derive the stellar 
masses of LBGs, LAEs and radio galaxies, and perform 
a comparison between the samples in overdense regions 
and in the (GOODS) field. 

3.1. UV-optical Color- Magnitude Diagrams in RG 
Fields at z ~ 4 — 5 

In Fig. [5]we show the 3.6/im versus z-[3.6] diagram for 
our samples of galaxies at z ~ 4.1 in the field TN1338. 
Symbols indicate the (/-dropouts (open circles), the LAEs 
(stars), and the radio galaxy (large filled circle). Two 
LAEs were unfortunately confused beyond the possibil- 
ity of deblending. Of the remaining 10 LAEs, only three 
were detected at 3.6/im and we derive limits for the seven 
that were undetected. Also in Fig. [5j we show the distri- 
bution of z — 4 £?4 35 -band (B) dropouts from GOODS 
(grey circles) as selected by IStark et al.l (|2009f ) for com- 
parison. 

Next, in Fig. [6]we present the CMD of the U-dropouts 
(open circles) and the radio galaxy (filled circle) at 
z ~ 5.2 in the field TN0924. Unfortunately, three out 
of four of the LAEs in this field suffered from a large 
amount of confusion, and no attempt was made for pho- 
tometry. One LAE was found to be relatively isolated 
but not detected. Again, for comparison with the field 
we indicate a large number of U-dropouts from GOODS 
using grey circles. 

A number of basic observations from these diagrams 
can be made: 

(1) Our dropout samples show a sequence in the CMD 
with a total width in color of — 2 magnitude and a 
total range in 3.6/im of four magnitudes (~22-26 mag), 
peaking at 3.6/im — 25 and z-[3.6]~ 1. While the 
upper envelope in both CMDs is largely set by the 
limiting z-band magnitude by which the samples were 
selected (dashed lines in Figs. [5] and [5]), the lower 
envelope is physical in the sense that brighter objects 



(at 3.6/xm) are also redder. The same trend is seen for 
the GOODS reference samples (grey circles), indicating 
that it is a general feature of the high rcdshift population. 

(2) The LAEs in the TN1338 field, the only field with 
a sufficiently large sample available for attempting 
photometry by means of deblending (see §2), mostly 
show non-detections at 3.6/im. The three detecte d 
LAEs were L9, L16 and L22 (see lOverzier et al.l 120081 ) . 
L9 is also part of the LBG sample, while L16 and 
L22 are respectively too blue or too faint to make it 
into our (/-dropout LBG sample. These results were 
expected given that about half of the LAEs are fainter 
than Z85o=26.5 magnitude in the rest-frame UV and 
the rather linear correlation between the rest-UV and 
optical in the CMD. In any case we can conclude that 
we have not found any LAEs that are particularly bright 
or red compared to the LBGs, although we should note 
that the RG has an extremely luminous associated Lya 
halo typical of radio galaxies that would formally put it 
into a Lya-selected sample as well. 

(3) Both radio galaxies have extraordinarily bright 
rest-frame optical magnitudes, at least one magnitude 
brighter than the brightest LBG and 2-4 magnitudes 
brighter at 3.6/im compared to the bulk of the LBG pop- 
ulations at z ~ 4 and z ~ 5. They also have relatively 
red UV-optical colors but perhaps not redder than the 
extrapolated color of the sequence formed by LBGs at 
lower magnitudes. 

3.2. Mass Segregation of LAEs, LBGs and Radio 
Galaxies 

As explained in i j2.7.2l we use a limited set of spectral 
energy distributions as simple baseline models for inter- 
preting the CMDs shown in Figs. [5] and [6] In these 
figures we have indicated the expected location of star- 
forming galaxies with exponentially declining star forma- 
tion histories of r ~ 100 Myr, dust attenuations in terms 
of E(B - V) of 0.0 (blue lines), 0.15 (green lines) and 0.3 
mag (red lines) modeled using t he recipes for starburst 
galaxies given in ICalzettil (|2001h . and plotted at three 
ages of 100, 200 and 300 Myr (indicated by the circle, 
triangle and square, respectively). The tracks were nor- 
malized to a total stellar mass of 10 10 Mq. These models 
span an appropriate range of par ameter spac e com pared 
to th e detailed investigations of IStark et al.l fe.g. 120071 . 
l2009l and references therein) and Bouwc ns et alj (|2009h . 
The right and top axes of Figs. [5] and [6] show the atten- 
uation and stellar mass obtained for a model having a 
fixed age of 200 Myr. 

3.2.1. RGs 

At z = 4.1 and z — 5.2, the 3.6/im band probes the 
rest-optical at ~7000 and 6000A, respectively. At these 
waveleng ths, the light is believed to be large ly stellar 
m origin ()Zirm et alj|2003t ISevmour et al.ll2007D . despite 
the many significant, non-stellar components of radio 
galaxies across the entire spectrum. In both fields, we 
find that the RG is the most massive object, with an 
inferred mass of ~ 10 11 Mq. Without more constraints 
on the SED we judge that this estimate can be higher or 
lower by —25% depending on the exact age and atten- 
uation. TN0924 at z = 5.2 is over a magnitude redder 
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Fig. 5. — The 3.6/^m versus z— [3.6] color-magnitude diagram for galaxies in the field TN1338. Symbols indicate the radio galaxy (large 
filled circle), LBGs (open circles) and LAEs (filled stars). The 3.6^tm flux of the radio galaxy has been corrected for Ha emission line 
contamination (see !|2.6> . The S435 -dropouts from [Stark ct al. (200j|) are shown for comparison (small grey circles). The approximate 
detection limits for the TN1338 field and the GOODS control field are ind icated by the black solid and grey dotted lines, respectively. The 
colored solid lines indicate our BC03 baseline SED models described in t!2 . 7. 2 1 i.e. exponentially declining SFHs (t = 100 Myr) having 
attenuations of E(B — V) = 0.0 (blue), 0.15 mag (green) and 0.3 mag (red), each plotted at ages of 100 (circle), 200 (triangle) and 300 
Myr (square). The models were normalized to a total stellar mass of 10 10 Mq. The right and top axes, respectively, show the attenuation 
and stellar mass obtained for a model having a fixed age of 200 Myr (for M„ we fix E(B — V) to 0.15). The red dotted line (marked 
'dust') indicates the minimu m color change as a f unction of 3.6/xm magnitude that is expected purely due to dust based on the M\jy vs. 
E{B — V) relation found by Bouwens ct al. (2009) (and further assuming r = 100, t = 200 Myr in order to predict the 2-[3.6] color). The 
blue dotted line (marked 'dust+age') indicates a toy model in which we let both the dust and age increase in steps of AE(B — V) = 0.05 
and At = 50 Myr with decreasing 3.6^tm magnitude starting from E(B — V) = 0.0 and t = 200 Myr (each step is indicated by a small blue 
square). The latter model ('dust+age') is at least qualitatively in better agreement with the relatively steep change in color as a function 
of the 3.6/mi magnitude compared to the 'dust only' model. 



than TN1338 at z = 4.1. If the attenuations are similar, 
it implies that the higher rcdshift RG is about twice as 
old as the lower one (~600 vs. 300 Myr). In any case, 
it is interesting that both RGs are much more massive 
than any of the LBGs in their surrounding fields. 

3.2.2. LBGs 

We now compare the dropouts with the model tracks 
plotted in Figs. [5]and[6l In both RG fields we find that 
the range in 3.6/xm magnitude of the g-dropout and V- 
dropout LBGs implies a distribution in stellar mass of 
about one order of magnitude, from a few times 10 9 to 
a few times 10 10 Mq. The CMD slope can be explained 
entirely by a correlation between dust and stellar mass, 
finding ages of ~ 100 — 200 Myr and dust increasing 
from ~ 0.0 to ~0.15 mag between the low and high mass 



end. Alternatively, interpreting the CMD slope as an 
age effect, the low mass end objects are consistent with 
an age of ^70-200 Myr and small attenuation, while the 
typical age at the high mass end would be ~ 2.5 x higher 
(~ 200 — 500 Myr) for a similar attenuation. We will 
discuss the age-dust (and metallicity) degeneracy in more 
detail in Sect. 13.31 below. 

3.2.3. LAEs 

Due to the non-detections and overall small number of 
LAEs in the TN0924 field, we can only compare LAEs 
and LBGs in the TN1338 field. Most of LAEs were unde- 
tected, suggesting masses of < 10 9 M©, and we note that 
many of the faintest LBGs were undetected as well. The 
brighter, detected LAEs have colors that are consistent 
with those of non-Lya LBGs at similar optical magni- 
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Stanwav et all 120051 : lHathi et all l2008t iBouwens et"all 
2009), in some contrast with earlier results that indi- 



FlG. 6. — The 3.6/xm versus z-[3.6] color-magnitude diagram for 
galaxies in the field TN0924. Symbols indicate the radio galaxy 
(larg e filled circle) and LBGs (open circles). The Vgoe-dropouts 
from lStark et ail l|200Sft are shown for comparison (small grey cir- 
cles). The approximate detection limits for the TN0924 field and 
the GOODS control field are indicated by the black solid and 
grey dotted lines, respectively. The baseline SED model tracks 
are shown for z ~ 5 but are otherwise identical to the models in 
Fig. [5] (see caption of Fig. [5] for details). The right and top axes, 
respectively, show the attenuation and stellar mass obtained for a 
model having a fixed age of 200 Myr (for M, we fix E(B — V) to 
0.15). The red dotted line (marked 'dust') indicates the minimum 
color change as a function of 3.6/im magnitude that is expected 
pu rely due to dust based on the Mjjv vs. E{B — V) relation found 
by IBouw ens et al. (2009) for z ~ 5. The blue dotted line (marked 
'dust+age') indicates a toy model in which we let both the dust 
and age increase in steps of AE(B — V) = 0.05 and At = 50 Myr 
with decreasing 3.6/xm magnitude starting from E(B — V) = 0.0 
and t = 200 Myr (each step is indicated by a small blue square). 

tudes. The tendency of LAEs at z ~ 4 to be relatively 
faint and blue in the rest-frame UV/optical is consis- 
tent with their young age as found in other surveys (e.g. 
Shaplev et al.|[2003tlGawiser et a l. 2006; Pcnt ericci et al.l 
20071 120091: lOuchi et alj 120081) ari d predicted by some 
simulations (jNagamine et al-lfeOOSfl . 



3.3. Interpreting the relation between and z-[3.6] 

color 

In Figs. [5] and [6] we have observed a strong trend of 
redder z-[3.6] colors with increasing stellar mass. This 
trend is qualitatively similar to the "blue envelope" o f 
LBGs at z ~ 3 - 4 observed by iPapovich et"aT1 (|2004D 
and must be related to an increase in the mean age, 
dust or metallicity (or a combination thereof) with in- 
creasing rest-frame optical luminosity. While we can- 
not solve these degeneracies on an object-by-object ba- 
sis without deep multi-wavelength photometry, we can 
use recent results from literature to explain the gen- 
eral trends observed. Deep studies of dropouts covering 
a wide range in both redshift and UV luminosity have 
now clearly demonstrated that star-forming galaxies be- 
come bluer in the UV continuum with increasing redshift 
and with decreasing luminosity (|Papovich et al.l 12004 



cated no clear correlation between t he slope of the UV 
continuum and UV luminosity (e.g. iMeurer et alJH999t 
lAdelberger fc Ste idel 2000: lOuchi etaiT 2004). As shown 
bv IBouwens et all (|2009r) . these results can be fully ex- 
plained by the relatively limited range in redshift and/or 
UV luminosity probed by earlier surveys. Using the 
deepes t data available from t he Hubble Ultra Deep Field 
(UDF; iBeckwith et al.l I2006D . these authors show that 
there is a significant correlation between the UV slope 
and the UV luminosity for both B- and U-dropouts. 
They also investigate the relative effects that, in prin- 
ciple, metallicity, age and dust each could have on the 
observed UV slope. They find that both dust and age 
have the potential of reddening the UV continuum sig- 
nificantly, with the effect of dust being the strongest for 
realistic distributions of age and dust. The influence of 
metallicity is found to be only small. 

Could the trends that we found in Figs. [5] and [5] 
be explained due to systematic changes in dust and 
age with in c reasin g stellar mass (or 3.6/im luminosity)? 
IStark et al.l (|2009t ) have shown that the median ages of 
dropouts at 4 < z < 6 (the same samples that we have 
used for our comparison) change only little with either 
redshift or UV luminosity (less than a factor of 2). How- 
ever, they do find a strong correlation between stellar 
mass and SFR as inferred from the UV continuum lumi- 
nosity (with a net specific SFR, SFR/M*, that is rela- 
tively constant as a function of stellar mass). The corre- 
lation between SFR and extinction by dust is well-known, 
as galaxies with higher SFRs tend to have more dust (e.g. 
Wang fc Heckman 1996; Ma rtin et all2005l : iReddv et al.l 



20061) . IBouwens et all (|2009l ) have now firmly established 
this correlation for LBGs by studying the average UV 
color (a measure of reddening) as a function of the abso- 
lute UV magn itude (a measure of th e unextincted SFR). 
As argued bv IBouwens et al.l (|2009h . we thus expect to 
see a strong correlation between stellar mass and redden- 
ing due to dust. In Figs. [5] and [5] we have indicated the 
expected amount of reddening purely due to dust based 
on the empirical UDF results. In both figures, the red 
dotted line marked 'dust' indicates the color change (as 
a function of 3.6/xm magnitude) expecte d based on the 
relatio n between M\jy vs. E(B — V) from lBouwens et al.l 
(2009, see their Fig. 3). Using this relation, the extinc- 
tion varies from E(B — V) «0.0 to psO.2 mag from the 
faint to the bright end. While the trend predicted by 
our UV-to-(z-[3.6]) color extrapolation shows the cor- 
rect general behavior, the reddening in z-[3.6] expected 
seems to underpredict the amount of reddening observed. 
This may indicate that, besides the increased reddening 
due to dust, additional reddening is required due to age 
as a function of stellar mass. 

To illustrate this further we construct a simple toy 
model. We let both dust and age increase towards 
brighter 3.6^m magnitudes. The E(B — V) is increased 
in steps of 0.05 and the age in steps of 50 Myr. The 
resulting color trend for this toy model is shown by the 
blue dotted line marked 'dust+age' in Figs. [5] and [6l 
The starting point (marked by the small blue square on 
the faint end of the blue dotted line) was chosen to have 
E(B -V) = 0.0 and an age of 200 Myr. The 'dust+age' 
model is, at least qualitatively, in better agreement with 
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the relatively steep change in color as a function of the 
3.6/im magnitude compared to the 'dust'-only model. 

In principle, it should be possible to find the required 
evolution in dust and age that best fits the general trends 
observed in the CMDs. However, Figs. [5] and [6] show 
that the dynamic range in z— [3.6] color in the TN1338, 
TN0924 and GOODS fields is too small to carry out such 
an exercise (i.e., a fit to the data can not be made due 
to relatively red sources being missed at the faint end). 
Deeper data in the rest-frame UV would greatly alleviate 
this problem, and it should be possible to perform such 
a study in the much deeper UDF. 

3.4. Comparison of galaxies in overdense regions versus 
galaxies in average environments from GOODS 

It has previously been shown that the number den- 
sity of dropouts and LAEs in the TN1338 and TN0924 
fields is larger relative to random fields, indicating 
that these fields may host massive gal a xy clusters-in- 
formation ([Venemans et al.l I2002L l2004t lOverzier et al.l 
120061 12008D . One of the goals of the current paper was 
to test whether besides this enhancement in number den- 
sity there also exist differences in the physical properties 
of galaxies in overdense regions compared to the field. 
The large samples of B- and ^-drop outs extracted from 
the ~ 320 arcmin 2 GOODS field by IStark et alj (gflOg ) 
allow us to make such a comparison. Comparing the 
color-magnitude distribution in our fields with the re- 
sults found for GOODS (grey circles in Figs. [5] and [6]) 
indicates a very similar trend of redder colors for brighter 
objects. This suggests that, at least at first sight, in both 
types of environment the dropout populations as a whole 
possess very similar stellar populations. 

Unfortunately, there are several factors that compli- 
cate making further quantitative statements about the 
differences or similarities between overdense fields and 
GOODS. First, the size of the individual error bars in 
the CMD diagrams is substantial, thereby limiting the 
usefulness of object-by-object comparisons except for the 
brightest objects (as we will see below). Second, because 
we only have limited filters and miss deep NIR data, an 
analysis based on multi-band SED fitting is not yet pos- 
sible for our fields. Third, although many of the dropout 
galaxies are believed to be associated with the protoclus- 
ters in the vicinity of the radio galaxies, we must note 
that the dropouts in TN1338/TN0924 are not spectro- 
scopically confirmed, except for the subset that are Lya- 
bright. This means that our dropout samples will un- 
doubtedly include galaxies in the (near) fore- and back- 
ground of the protoclusters that could offset or bias any 
trends observed in Figs. [5] and [6] 

These caveats aside, however, we can at least ask the 
question what differences with GOODS we would in 
principle be able to detect based on the quality of our 
data presented. The first possibility we consider is that 
if galaxies in the overdense regions were systematically 
bluer or redder compared to the field, the protoclustcr 
dropouts may be expected to form a lower or upper en- 
velope with respect to the well-established CMD slope 
formed by the GOODS dropouts at z ~ 4 — 5. We esti- 
mate that a systematic difference in ages by a factor of 
> 2 or a difference in attenuation of > 0.15 mag would 
have been detectable given the size of the error bars in 
Fig. [5] However, no such differences are suggested by 
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[3.6] 

Fig. 7. — 3.6/^m versus z-[3.6] color-magnitude diagram for the 
z ~ 6 «775-dropouts selected from the fields SDSS0836 (QSO at 
z = 5.82) and SDSS1030 (QSO at z = 6.28). Symbols indicate 
the quasars (large filled circles) and i775-d ropouts (open circles). 
The 4775-dropouts from Stark et al. (2009) are shown for compar- 
ison (small grey circles). The approximate detection limits for the 
QSO fields and the GOODS control field are indicated by the black 
solid and grey dotted lines, respectively. Model tracks have been 
indicated for z ~ 6 but are otherwise identical to the ones shown 
in Figs. [5]and Fig. [6] (see the caption of Fig. [5]for details). The 
QSOs have been indicated for reference only, as their magnitude 
and colors are completely dominated by AGN light instead of stel- 
lar light. 

the data. 

Another possibility we consider concerns the hypoth- 
esis that the overdense regions contain more massive 
galaxies compared to the general field. In this scenario, 
we would expect an extension towards brighter 3.6/im 
magnitudes (larger stellar masses) relative to GOODS. 
Because every object in the overdense region would be- 
come a little brighter, this means that at the faint end 
of the luminosity/mass function more objects in the 
overdense region would pass the detection threshold, 
while the number for the "field" would remain the same. 
Therefore, one would not immediately expect to find any 
differences between the average properties of galaxies at 
the faint end. At the bright end, however, one would 
expect an excess of massive sources relative to the field. 
We estimate that a boost in masses by a factor of a few 
should be detectable in the CMD diagrams. In this re- 
spect, it is interesting that in both fields, the RGs are 
brighter than any dropout from the GOODS compari- 
son sample, even though the latter survey covers an area 
that is over 25 times larger. We have also analyzed the 
relative numbers of dropouts in the fainter magnitude 
bins, finding no evidence for a significant brightening of 
dropouts in the overdense regions. From this it follows 
that compared to GOODS the protocluster fields likely 
do not contain an excessive number of bright/massive 
(~ 10 11 Mq) LBGs other than the radio galaxies them- 
selves. We will further discuss these results in 2J 
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3.5. UV-optical Photometry of i- dropouts in Two QSO 
Fields at z ~ 6 

For completeness, we will use our ACS and IRAC data 
on the two z ~ 6 QSO fields in order to perform a similar 
analysis for the «-dropout candidates in the (projected) 
vicinity of the luminous quasars. We have tabulated our 
measurements in Table 4 and plot the CMD in Fig. [71 
We also plot the QSOs in the CMD but stress that their 
emission is strictly dominated by non-stellar processes 
and should not be compared with any other populations 
or models shown in this paper. We compare our re- 
sults again with a small sample of similarly selected i- 
dropouts from GOODS (grey circles), and with the same 
model tracks as plotted in Figs. [5] and [6] The colors of 
the four i-dropouts in the QSO fields that are detected 
in both bands are clearly inconsistent with very young, 
dust-free models (e.g. the blue track in Fig. [Jj, while 
the masses are on the order of 10 10 M Q relatively inde- 
pendent of dust and age. The faintest objects that can 
be probed in GOODS are about ten times less massive 
due to the much higher sensitivity in the [3.6] channel 
(jYan et al.ll2005t IStark et al.ll2007f ). Comparing the two 
distributions further suggests an excess of sources hav- 
ing 3.6/im<24.5 magnitude in the QSO fields not seen in 
GOODS. Unfortunately, the small number statistics, the 
uncertainties in the redshifts and the photometric errors 
do not allow us to say at present whether this trend is 
real or not. 

4. SUMMARY AND DISCUSSION 

4.1. Summary 

Summarizing our results from the previous section, at 
z ~ 4.1, the masses and stellar populations of RGs and 
LBGs are found to lie in a similar range as those at 
z ~ 5.2, with the RGs being systematically the most mas- 
sive objects (~ 10 11 Mq). The stellar mass of the LBGs 
ranges from a few times 10 9 to a few times 10 10 M Q , and 
there is a clear trend of redder colors with increasing 
mass. The correlation can be explained if more massive 
galaxies have more dust (because their SFR is higher) 
and are somewhat older compared to less massive galax- 
ies. Few LAEs were detected, implying they are generally 
of low mass. The RGs, LBGs and perhaps LAEs thus ap- 
pear to be on different evolutionary tracks with the RGs 
having experienced either a more prolonged or more effi- 
cient phase of star formation history compared to LBGs, 
and LAEs being on average in earlier stages than the typ- 
ical UV-selected LBGs in these fields (but perhaps not 
different compared to LBGs of the same UV luminosity!). 
We investigated possible differences between the proper- 
ties of galaxies in overdense regions and those of dropouts 
from the much larger GOODS survey based on the distri- 
bution of objects in the rest-frame UV-optical CMD dia- 
gram. We found no significant, systematic differences be- 
tween the overall distributions in z-[3.6] color and 3.6/im 
magnitude between dropout galaxies in the different re- 
gions. We did however find that in both fields the stellar 
masses of the radio galaxies as inferred from the 3.6/im 
flux are highly excessive compared to the general field 
expectations from GOODS. The latter survey virtually 
contains no dropout galaxies as bright (at 3.6/im) as the 
RGs, even though the area surveyed by GOODS is about 
25 x larger. Their (inferred) stellar masses of ~ 10 11 M© 



are also about an order of magnitude larger than the 
brightest LBGs in the protocluster fields. 

4.2. The most massive galaxies at z ~ 4 — 5 

Even though we have shown that none of the objects 
in GOODS are as bright as the RGs, IStark et all (|2009l ) 
find a non-negligible number of objects having compara- 
ble masses based on the detailed SED fits performed for 
GOODS. Thus, for the first time, we have a quantitative 
number density of massive galaxies from the field to be 
compared with RGs and protocluster regions. In order 
to quantify how rare galaxies with stellar masses as high 
as observed for the RGs are, we use the total number 
densities of (UV-selected) galaxies having st ellar masses 
> 10 11 M Q as derived by IStark et all (|2009D . They find 
number densities of N(> M) ~ 6 x 10~ 6 - 4 x KT 5 
Mpc" 3 at z ~ 4 and N(> M) ~ 1 X KT 6 - 2 x KT 5 
Mpc -3 at z ~ 5 (ranges are la), and we note that 
the upper values on these number densities are perhaps 
overestimated as non-spectroscopically confirmed sources 
having z— [3.6] > 2 may have some con tamination from 
dusty star-formin g systems at z ~ 2 (|Yan et al.l 120041 : 
IStark et al.ll2009h . 

The co- moving volumes of the TN1338 and TN092 4 
fields are « 2 x 10 4 Mpc 3 (|Overzier et al.ll2006l 120081) . 
implying that we would naively expect as many as, re- 
spectively, 0.1-0.8 (0.1-0.6 if we assume that the 24/im- 
detected objects in GOODS are foreground contamina- 
tion) and 0.02-0.4 highly massive galaxies in a random 
single ACS pointing. Thus, it appears that the pres- 
ence of a single ~ 10 11 M (radio) galaxy in the pro- 
tocluster fields is marginally consistent with the (la) 
upper limits on the cosmic average number density at 
z ~ 4 while it is a positive deviation of several sigma 
at z ~ 5. The GOODS field estimates are about two 
times higher than the values derived by iMcLure et al.l 
(2008) based on a much larger survey of 0.6 deg 2 , indicat- 
ing that the GOODS estimates are probably not missing 
a significant fraction of massive (UV-selected) objects. 
Furthermore, to date no objects have been found with 
masses in excess of a few times 10 11 M Q . Our results 
thus confirm earlier studies that suggest that RGs are 
amo ng the most massive galaxies in the early Universe 
(e.g. iPentericci et al.lll997t ISevmour et al J 120071 and see 
Fig. 5 in iRocca-Volmerange et all (2004)). The low av- 
erage number density of such objects further implies that 
the RGs in the TN1338 and TN0924 fields have experi- 
enced a rapid growth through accretion of gas or merging, 
presumably as a consequence of their overdense environ- 
ments. This is qualitatively consistent with some models 
that specifically require an extra bias factor due to merg- 
ers in order to explai n the strong clust ering of luminous 
quasars at z — 4 (e.g. IWvithe fc Loebll2009» . 

From the lack of any known brigh t radio sources among 
the massive GOODS galaxies (see iMainieri et aTl 120081 
for details on a radio follow-up in the GOODS-CDFS), 
we can deduce that if a powerful (radio) AGN is a typ- 
ical phase of massive galaxies at z ~ 4 — 5, then it is 
unlikely to last longer than 30-50 Myr. Otherwise, ra- 
dio galaxies would be quite common in a field the size of 
GOODS given that the -B-dropout (U-dropout) samples 
span ~500 (^300) Myr in look-back time. This simple 
estimate based on the number counts is quite close to 
the expected radio source life-times of ~10 Myr derived 
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using other arguments (Blundel l fc Rawhngslll999T ). If, 
on the other hand, none of the massive GOODS galaxies 
become radio galaxies even for a brief period during the 
epoch at z ~ 4 — 5, it may mean that additional require- 
ments for producing them must be met. Perhaps the 
overdense environment as observed in the protocluster 
fields TN1338 and TN0924 facilitates the rapid accre- 
tion of cold gas th rough mergers or cosmological flows 
(|Keres et al.l I2009T) spawning both a massive starburst 
and powerful AGN activity of a supermassive black hole 
(see iZirm et alJl2005l l2009h . 

4.3. Relation to low redshift clusters 

iStark et all ([20091) compared the number density of 
> 10 11 M Q B— and V-dropout systems to the observed 
density of quiescent ~ 10 11 M© galaxies at z ~ 2.5, find- 
ing that the high redshift sample can account for ~ 20% 
of the massive population at lower redshift. In a simi- 
lar fashion, we can compare the number density of the 
most massive objects in GOODS at these redshifts to the 
number density of galaxy clusters in the local Universe, 
finding that they are roughl y equal for cluster s having 
L x ~ 1 - 5 x 10 43 erg s" 1 (jRosati et al1ll998l ). Thus, 
one may expect that the z = descendants of the most 
massive z ~ 4 galaxies in GOODS are fo und in typi- 
cal gr oup and cluster environments (see also lOuchi et al.l 
l200l . Since the radio galaxy fields studied in this paper 
are overdense with respect to GOODS, it will be inter- 
esting to perform a systematic study of galaxies in the 
direct environmen t of the massive GOODS galaxies (e.g. 
iDaddi et alJl200flh . 

We have argued that the overdense regions correspond 
to sites of massive cluster formation. If so, we expect 
that the protocluster region will start to develop a clus- 
ter red sequence and a virialized intracluster medium at 
some redshift z < 2. Likely, the stars formed in the 
protocluster galaxies observed at z ~ 4 — 5 will evolve 
and become incorporated into the population of massive 
early-type galaxies on the red sequence. We can use our 
results of Sj3]to estimate what fraction of this low redshift 
cluster red sequence mass is represented by all protoclus- 
ter galaxies detected. We use the conversion between 
3.6/im magnitude and M* derived from the simulations 
in §2.7.31 and calculate the total mass of all dropouts in 
the TN1338 field at z = 4. We find a total stellar mass 
of ~ 5 x 10 11 Mq. About 25% of this mass is due to the 
large contribution from the radio galaxy. 

We will now compare our estimate to the total stellar 
mass on the red sequence of the cluster RDCS1252.9- 
2927 (C1125 2), a massive X-ray l uminous cluster at 
z = 1.24 (iBlakeslee et al.1 l2003at iRosati et all 12004 
iDemarco et alll2007j ). Although we do not necessarily 
expect that TN1338 will evolve into a similar cluster as 
C11252 (see discussion below), we choose this particular 
cluster for our demonstration for a number of reasons. 
First, the analysis of C11252 was performed by members 
of our team, giving a good understanding of the data and 
its interpretation. More importantly, it is one of the best 
studied high redshift clusters, and its high redshift en- 
sures a minimal gap in cosmic time with respect to z = 4 
(about 3.5 Gyr). This is important because below we 
will perform a simple extrapolation of the stellar masses 
measured at z ~ 1 to those predicted at z = 4. 

The analysis of C11252 performed by iRettura et al.l 



(2008) found that 80% of the ETGs have population 
ages of 3.5 ± 1 Gyr, consistent with a mean formation 
redshift of z ~ 4. If we limit ourselves to ETGs for 
which the mean star formation weighted ages are > 3.5 
Gyr, we find 10 massive (M* > 10 11 Mq) early-type 
galaxies (ETGs) representing a total "red sequence 
mass" of 3.5 x 10 12 Mq (about half of which is due 
to the three most massive ETGs). If we now assume 
that all the dropouts detected in the TN1338 field are 
part of the protocluster, its total stellar mass amounts 
to about 14% (at most) of the total C11252 ETG mass 
presumed to have formed at z > 4. This percentage 
(14%) implies that we are currently missing most of the 
mass that we naively expected to have existed already 
at z ~ 4. Where could this mass be hiding? We suggest 
four possible explanations for this apparent discrepancy: 

(1) When evolved to z ~ 1, TN1338 could be a 
very different cluster than C11252. On one hand, large 
uncertainties in the physical sizes and the magnitude 
of the (dark matter) overdensities make it currently 
very difficult to predict the total mass and virialization 
redshift of the descendants of protoclusters. Estimates 
indicate that their total masses range from 10 14 to 
10 15 Mq, with virializ a tion redshifts of < z < 1 
(e.g. ISteidel et all 120051 : lOverzier et all I2008D . On the 
other hand, C11252 is an exceptionally massive, X-ray 
luminous cluster indicating virialization that was largely 
complete at z = 1.24. The relatively low total stellar 
mass found for TN1338 could therefore indicate that it 
will evolve into a more typical and lower mass cluster, 
containing fewer ETGs on the red sequence at z < 1 
compared to C11252. Interestingly, the existence of 
extremely massive clusters such as C11252 at z ~ 1.2 
may indicate that much larger protoclusters at z > 2 
remain to be discovered. 

(2) A large fraction of the ETG mass may be accreted 
from a much larger region that extends beyond the region 
pr obed by our ACS pointing. The Lya emitters found 
by IVenemans et al.l (|2007| ) extend signific antly beyond 
the ar ea covered by our HST data, and llntema et al.l 
(2006) found evidence that TN1338 is part of a larger 
structure of LBGs at a co-moving scale of ~10 Mpc. 

(3) The stars may be contained in numerous sub- 
units that are too faint to be detected by our current 
survey. Extending the luminosity function of z ~ 4 
-B-dropouts down to significantly below the detection 
limit suggests that the contribution from faint objects 
to the total star form ation rate amounts to <50% (e.g. 
iBouwens et all 120071 ). We may therefore increase our 
total mass estimate in the overdense region by no more 
than about a factor of 2, increasing the ETG progenitor 
mass accounted for at z ~ 4 to ^30%. If this very faint 
high-redshift population provides indeed an important 
contribution to the mass of ETGs, they will need to 
merge to form the ETG population within a few Gyr. 

(4) A significant fraction of the stellar mass may be 
contained in systems that are either very dusty or have 
little on-going star formation and are hence missed by 
our UV-selection. 
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The answers to questions (1) and (2) are currently dif- 
ficult to address based on the limited observational data 
and the small number of protoclusters known. These 
issues will be the subject of a future paper in a se- 
ries of papers aimed at comparing the observational 
properties of statistical samples of (proto-)clusters ex- 
tract ed from cosmological simulations of cluster forma- 
tion ([Overzier et alj [20091 and in prep.). The possible 
incompleteness due to red galaxies (question 4) will be 
discussed in the next section. 

4.4. Have we missed a significant population of 
quiescent or obscured galaxies? 

Because UV-selected samples may miss galaxies that 
are either quiescent or obscured, we may expect that the 
estimates for the number density of massive galaxies both 
in GOODS and in the protocluster fields provide a lower 
limit to the actual number densities. At z > 2 there 
is a substantial contribution to the stellar mass func- 
tion from the populations of so-called distant red galaxies 
(DRGs), which consist of both obscured, star- forming ob- 
jects and red, quiescent objects that tend to be substan- 
tially more massive (and older) compared to UV-selected 
LBG s (e.g . Franx et al. 2003; Brammcr & v an Dokkuml 
[20071: IDunlop et all [20071 : iWiklind et aT][2008h . If such 
objects are abundant in protocluster regions, we may 
largely resolve the discrepancy we found above between 
the total stellar mass observed and that expected based 
onz ~ 1 ETGs in C11252. Will the overdensities in 
LBGs and LAEs observed in TN1338 and TN0924 be 
accompanied by simil ar overdensities in relativ e ly red 
or massive galaxies? iBrammer fc van Dokkuml (2007) 
have shown that the UV-selection at z ~ 2 — 3 may 
miss as much as ~ 33% of galaxies having very red 
UV continuum slopes that are selected on the basis of 
a strong Balmer Break. However, at a higher redshift of 
z ~ 3 — 4.5 almost all (> 90%) galaxies in their Balmer 



Brea k sample are als o selec ted by means of a UV selec- 
tion. iBouwens et al.l (|2009h show that the UV-dropout 
selections in GOODS and the UDF at z > 2 become sub- 
stantially more complete with increasing redshift. This 
is because the UV-continuum slope of the galaxy pop- 
ulation becomes bluer as redshift increases. At z > 5 
the distribution of UV slopes of the entire dropout pop- 
ulation lies within the dropout selection function, and 
even at z ~ 4 the selection is believed to be largely 
compl ete co nsistent with, e .g.. IBrammer fc van Do kkum 
(|2007ft and IDunlop et at] ([2007J). Finally, we consider 
the population of ultra- luminous IR galaxies (ULIRGs). 
The high redshift ULIRG population is a population 
of heavily obscured starbursts detected in the sub-mm 
and corresponds to rapidly forming, massive galaxies at 
z ~ 2 — 4. However, analogous to the population of 
Balmer break objects, the contribution of ULIRGs to the 
star formation rate density declines with redshift from 
~30% at z ~ 2.5 t o ~ 10% at z ~ 4 (ICapak et all 120081 : 
iDaddi et al.1 120091 : IBouwens et al.ll2009f l. Therefore, we 
currently do not have very strong reasons for believing 
that a large population of either the red Balmer break 
objects or heavily enshrouded starbursts is being missed 
in our TN1338 protocluster at z = 4. Follow-up studies 
are currently being undertaken (C. De Breuck, private 
communication) based on deep near-infrared data specif- 
ically targeting galaxies having prominent Balmer breaks 
between the H and K bands (K and 3.6/im) for z ~ 4 
(z ~ 5 — 6) that could resolve whether protoclusters have 
an overabundance of such red objects compared to the 
general field. 
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TABLE 1 
Overview of Targeted Fields 
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Redshift 


«J2000 


<5j2000 


Type 


Galaxy Samples 
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4.11 
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-19:42:30.5 
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g-dropouts, Lya emitters 


TN J0924-2201 


5.20 


09:24:19.89 


-22:01:41.3 


Radio Galaxy 


V-dropouts, Lya emitters 


SDSS J0836+0054 


5.82 


08:36:43.87 


+00:54:53.2 


Radio-Loud Quasar 


i-dropouts 


SDSS J1030+0524 


6.28 


10:30:27.09 


+05:24:55.0 


Radio-Quiet Quasar 


i-dropouts 
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TABLE 2 

Photometric Properties of 3475-Dropouts in the Field of TN J1338-1942. 
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a IDs refer to Overzier et al. (2008). b Kron magnitude. Spitzer IRAC 3.6 fim and 4.5 /xm magnitude 
from GALFIT. d The Spitzer confusion classes: (1) isolated and detected; (2) isolated but undetected; 
(3) confused, but GALFIT may help; (4) heavily blended. Lya Emitters. 1 The 3.6 /im measurement 
most probably includes a large contribution from the Ha emission line halo of the radio galaxy. Applying 
a correction for Ha yields a 3.6 fim continuum magnitude of 21.2 ± 0.1 (see §3.1 for details). This is the 
value plotted in Figs. [5] and [6] 
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TABLE 3 

Photometric Properties of Spectroscopically Confirmed Lyo Emitters in the Field of TN 

J1338-1942. 



ID a 


«J2000 


Sj2000 


*775 - ^850 


Z 850 


3.6/im c 


A.5fj,m c 


Class™ 
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13:38:22.46 


-19:43:33.68 
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0.00 ±0.12 


25.81 ±0.15 


> 25.2 
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2 



a IDs refer to Overzier et al. (2008). b Kron magnitude. Spitzer IRAC 3.6 jj,m and 4.5 /jm mag- 
nitude from GALFIT. d The Spitzer confusion classes: (1) isolated and detected; (2) isolated but 
undetected; (3) confused, but GALFIT may help; (4) heavily blended. 



TABLE 4 

Photometric Properties of V606-Dropouts in the Field of TN J0924-2201. 
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a IDs refer to Overzier et al. (2006) . b Kron magnitude. Spitzer IRAC 3.6 
/im magnitude from GALFIT. d The Spitzer confusion classes: (1) isolated 
and detected; (2) isolated but undetected; (3) confused, but GALFIT may 
help; (4) heavily blended. c Lya Emitters. 



Stellar masses at z — 4 — 6 



TABLE 5 

Photometric Properties of «775-Dropouts in the Field of SDSS0836+0054 and 

SDSSJ1030+0524. 
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a IDs refer to Zheng et al. (2006) and Kim et al. (2008) . b Kron magnitude. Spitzer 
IRAC 3.6 /j,m from GALFIT. d The Spitzer confusion classes: (1) isolated and detected; 
(2) isolated but undetected; (3) confused, but GALFIT may help; (4) heavily blended. 



